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^ ■ Abstract 

H ; 

We study how to probe the anomalous CP- violating couplings of top quark with 7, Z 
and g in the tt threshold region at future e^e~ colliders. These couplings contribute to 
the difference of the t and i polarization vectors 6P and to the CP-odd spin correlation 
tensor 6Q,ij. We find that typical sizes of dP and 5Qij are 5-20% times the couplings 
{dfy, diz, dtg) in the threshold region. Experimentally 5P can be measured efficiently using 
the CP-odd combination of the £^ momenta or of the £^ directions. We have similar 
sensitivities to both the real and imaginary parts of the couplings independently using 
the two components of 6P. Taking advantage of different dependences of 5P on the 
polarizations and on the cm. energy, we will be able to disentangle the effects of the three 
couplings dt^ , dtz , dtg in the tt threshold region. We give rough estimates of sensitivities 
to the anomalous couplings expected at future e+e~ colliders. The sensitivities to dtj 
and diz are comparable to those attainable in the open-top region at e^e~ colliders. 
The sensitivity to dtg is worse than that expected at a hadron collider but exceeds the 
sensitivity in the open-top region at e'^e~ colliders. 
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1 Introduction 



Among all the fermions included in the Standard Model (SM), the top quark plays a very 
unique role. The mass of top quark is by far the largest and approximates the electroweak 
symmetry breaking scale. In fact the top quark is the heaviest of all the elementary particles 
discovered up to now. It means that in the SM Lagrangian the top quark mass term breaks the 
SU{2)i X U{1)y symmetry maximally. This fact suggests that the top quark couples strongly 
to the physics that breaks the electroweak symmetry. It is therefore important to investigate 
properties of top quark in detail, for the purpose of probing the symmetry breaking physics 
as well as to gain deeper understanding of the origin of the flavor structure. The standard 
procedures for investigating top quark properties are: measurements of fundamental quantities 
such as its mass and decay width; detailed examinations of various interactions of top quark 
to see if there are signs of new physics. Among them testing the CP-violating interactions 
of top quark is particularly interesting. This is because: (1) CP-violation in the top quark 
sector is extremely small within the SM. If any CP-violating effect is detected in the top 
sector in a near-future experiment, it immediately signals new physics. (2) There can be many 
sources of CP-violation in models that extend the SM, such as supersymmetric (SUSY) models, 
Leptoquark models (including P-parity violating SUSY models), multi-Higgs-doublet models. 
Extra-dimensions, etc. (3) In relatively wide class of models beyond the SM, CP-violation 
emerges especially sizably in the top quark sector. 

Predictions of certain models are as follows. In the SM, the lowest-order contributions to the 
electric-dipole-moment (EDM) of a quark come from three loop diagrams and are proportional 
to C^pas III]. Assuming that the results for u- and (i-quarks can also be applied to the top 
quark, one may estimate the top quark EDM as ~ 10~'^°ecm. One may also estimate the 
Z-EDM and chromo-EDM of top quark as ~ 10^^°ecm and ~ IQ^^^^ggCm, respectively, since 
there seems to be no reason that these two EDMs are much suppressed or enhanced. All 
these EDMs of top quark are quite small compared to those corresponding to the "(9(1)- 
couplings", e/mt ~ lO^^^ecm and Qs/mt ~ 10~^^(7sCm. On the other hand, the top quark 
EDMs are induced at one loop in many models, including multi-Higgs-multiplet models and 
SUSY models. In the two-Higgs-doublet models, a neutral Higgs can violate CP through 
the Yukawa interaction ^(a — 075)?/^ The size of the induced EDM is estimated asQ 

~ eCi^mf /(47r^m|) = 3 x lO^^^e cm (m^/lOO GeV)~^. The explicit calculations show somewhat 
smaller values, 10^^^-10^^°e cm, depending on ^/s and m^. In the Minimal SUSY Standard 
Model CP can be violated in the soft SUSY breaking sector |]6|, pl-pl . It was shown that the 
top quark EDM of ~ lO^^^ecm can be induced by gluino and chargino exchanges, assuming a 
universal gaugino mass and non-universal other soft-breaking parameters. 

Present experimental limits on the EDMs of top quark are not stringent [0. A limit on 
the chromo-EDM < lO'^^QgCm is obtained from the analyses using atot{pp — > tiX) and pt 
distributions of prompt photons produced in qg — > qj, etc. at Tevatron and from the analyses 
using Bt{B Xg'y) at CLEO. The limit on the EDM from the prompt photon distribution is 
similar: EDM < lO^^^e cm. 

* Here, one power of rrit is necessary to flip chirality. The extra two powers of TOj come from the Yukawa 
interaction. 
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There have been a number of sensitivity studies on the top quark EDMs expected at fu- 
ture hadron colhders and in the open-top region (y^ ^ 2mt) at future e+e~ and 77 coUid- 
ers. In hadron colhder studies 0,^,|T^,|T^-|16|, it is claimed that with the observablesf] made 
of elaborated combinations of momenta of charged leptons, 6-quarks, etc., experiments at 
^/s = 500 GeV and with an integrated luminosity 10fb~^ can probe the chromo-EDM down to 
a fewx 10~^''-10~^^(7s cm, or even to a fewx 10~^^ Qs cm by raising complexity of the observables. 
However, none of them performs detector simulations, which seem to be indispensable for a seri- 
ous sensitivity study.| Among several proposed CP-odd observables, lepton energy asymmetry 



= Ei+ — Ei- would be the simplest one [3,|T0|,|^. It is claimed that is sensitive to the 
imaginary part^ of the chromo-EDM down to ~ 10~^^gs cm assuming an acceptance efficiency 
e = 10%. Also studies of the EDM and Z-EDM of top quark in the open-top region at e~^e~ 
colliders are given in p|-p|,pn],pj|-|20[]. We take as a reference the results of which is based 



on simulation studies incorporating experimental conditions expected at a future e'^e~ linear 
collider. It is shown that, by using the mode tt — > bbWW — >• bbqq'ii' {i = e, /i), sensitivities to 
(the real and imaginary parts of) the EDM and Z-EDM are ~ 10~^^e cm at ^/s = 500 GeV, 
assuming an integrated luminosity of 10fb~^ and electron beam polarization of ±80%. Sen- 
sitivity to the top chromo-EDM in the open-top region at e^e~ colliders is studied in ||21|| ; 



they estimate a sensitivity ~ 10~^^(fe cm at y/s = 500 GeV, assuming an integrated luminosity 
50fb~^, an identification efficiency for top-pair production events ~ 100%, and E^^^ = 25 GeV. 
The last entry is a cut for the minimum gluon-jet energy, on which the sensitivity depends cru- 
cially. No detector simulation is performed in this study. The sensitivities of 77 colliders are 
studied in they are shown to be similar to those of e~^e~ colliders. 

Certainly it is desirable to probe the top quark anomalous interactions at highest possible 
energy where we have more resolving power, which motivated the above studies. On the other 
hand, it is known that studying various top quark properties in the tt threshold region at 
future e~^e~ colliders is promising and interesting; particularly the top quark mass will be 
determined to unmatched precision. A number of analyses elucidated physics potential of 



experiments in the tt threshold region [23-35]. Most of them, however, dealt only with the 



SM interactions. In this paper we extend these analyses and study how to probe anomalous 
CP-violating interactions of top quark in the ti threshold region. We note that there are some 
specific advantages in this region: 

• The polarization of top quark can be raised to close to 100% by adjusting longitudinal 
polarization of e~ beam ||36|,^. 

• Since top quarks are produced almost at rest, one can reconstruct the spin information of 
top quarks from distributions of their decay products without solving detailed kinematics 
131. 



The QCD interaction is enhanced in this region, so the cross section is sensitive to the 



t They include the so-called "optimal observables" |13| Q . 

^ For instance, it is important to study the effects of mis-assignment of jets to partons in event reconstruc- 
tions. 

§ Note that A^^ probes absorptive part of an amplitude A4, since it is CPT-odd. 
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top-gluon (tg) couplings. We can study anomalous tg couplings in a clean environment 
in comparison to liadron colliders. 

• There are less backgrounds from multiple W,Z productions compared to the open-top 
region. 

• In certain models (e.g. those in which a neutral Higgs is exchanged between t and t [§])) 
the induced top quark EDM and Z-EDM are enhanced near the tt threshold. 

Thus, for the sake of comparison with other kinematical regions, we would like to know sensi- 
tivities to CP-violation achievable in e'^e~' — > ti in the threshold region when these advantages 
are taken into account. 

In Sec. we present a qualitative picture of the effects of the anomalous CP-violating 
interactions in the threshold region. We derive the top quark vertices including the QCD 
enhancement in Sec. |^. The formulas for the polarization vectors and the spin correlation 
tensor of t and t are presented in Sec. followed by their numerical analyses in Sec. |^. 
Sec. P discusses the observables to be measured in experiments and gives rough estimates of 
sensitivities to the anomalous couplings. We summarize and conclude our analyses in Sec. |^. 
Some of the notations used in this paper are collected in the Appendix. 



2 Physical Picture 

Let us first review the time evolution of t and i, pair-created in e^e~ collision just below 
threshold, within the SM. They are created close to each other at a relative distance r ~ 1/mt 
and then spread apart non-relativistically. When their relative distance becomes of the order 
of the Bohr radius, r ~ (a^mt)"^, they start to form a Coulombic boundstate. When the 
relative distance becomes r ~ (mtFj)"^/^, where Tt is the decay width of top quark, either t or 
t decays via electroweak interaction, and accordingly the boundstate decays. Numerically these 
two scales have similar magnitudes, (asfrit)'^ ~ (mtFt)"^/^, and are much smaller than the 
hadronization scale ~ 1/Aqcd- Since gluons which have wavelengths much longer than the size 
of the tt system cannot couple to this color singlet system, the strong interaction participating 
in the formation of the boundstate is dictated by the perturbative domain of QCD. The spin 
and PC of the dominantly produced boundstate are J^*" = 1 . Inside this boundstate: t and 
t are in the S-wave state (L = 0); the spins of t and i are aligned to each other and pointing 
to e~ beam direction | tt ) or to e"*" beam direction \ H) or they are in a linear combination of 
the two states (S = 1). 

In this paper we consider anomalous CP-violating interactions of top quark with 7, Z, and 
g. In particular, we consider the lowest dimension effective operators which violate CP: 

jCcp-oM = -p^{Ua^^'i,t)d,A, - ^^{Ua^'^^,t)d,Z, - 3-^{tta^^-^,T'^t)d,Gl , 
2mt Init Irrit 

a^-^i[7^7^] , (1) 

where e = gw sin 9w and gz = gw / cos 9w- These represent the interactions of 7, Z, g with 
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Figure 1: The diagrams which contribute to the spin-dependent CP-violating potential between 
non-relativistic t and t. The vertex (gi represents the CP-odd interaction of top quark with gluon; cf. 
eq. (|l|) . An exchange of the Coulomb gluon g'~' gives the leading contribution of the CP-odd interaction 
to the potential. 



the EDM, Z-EDM, chromo-EDM of top quark, respectively]^ Each of these interactions has 
C = +1 and P = —1. We assume that generally the anomalous couplings dt^, dtz, dtg are 
complex where their imaginary parts may be induced from some absorptive processes beyond 
the SM. For a non-relativistic ti pair produced in e^e~ collision, the anomalous couplings of 
tt to 7 and Z reduce to 
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A ■ xK-^V)x* + ■ x\{-^'^)Xu (2) 

rrit rrit 

where Xt and Xt denote the two-component non-relativistic fields of t and t, respectively. The 
anomalous top-gluon coupling generates effectively a spin-dependent potential between t and i 

VcP-odd = ^{st- St) ■ VVcir) (3) 
rrit 

through the diagrams shown in Fig. |l]. Here, St and St denote the spins of non-relativistic t and 
t, respectively; Vc(r) = —Cpds/r is the Coulomb potential with the color factor Cp = 4/3. 
When dtg > 0, the potential Vcp-odd tends to align both chromo-EDMs in the direction of 
chromo-electric field, or, align St in the direction of r = — and St in the direction of — r. 

Let us consider effects of these anomalous interactions on the time evolution of the tt system. 
Assuming that the anomalous couplings dt^, dtz, dtg are small, we consider the effects which 
arise in linear perturbation in these couplings. CP-violation originating from the tj or tZ 
coupling occurs at the stage of the pair creation, i.e. when t and t are very close to each other. 
The generated boundstate has J^^ = , so t and t are in the P-wave (L = 1) and spin-0 
state I ti) — I it)- On the other hand, CP- violation originating from the tg coupling takes 
place after the boundstate formation when multiple gluons are exchanged between t and i, i.e. 
when t and t are separated at a distance of the Bohr radius. Therefore, first the boundstate 
is formed in J^'^ = 1 (L = and S = 1) state and after interacting via the potential 
Vcp-odd it turns into J^'^' = 1+^ (L = 1 and S = 0) state. Since we are interested in the 
dependences of observables on the couplings dt^, dtz, dtg up to hnear terms, we are interested 
in the interference of the leading SM amplitude and the amplitude including these couplings. 
The strong phases of these amphtudes that arise from QCD binding effects can be calculated 
reliably using perturbative QCD. 

^ The magnitudes of these EDMs are given by edtj/nit, gzdtz/iTit, gsdtg/mt, respectively. 
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Which CP-odd observables are sensitive to the above CP-violating couphngs? For the pro- 
cess e^e~ tt, we may conceive of following expectation values of combinations of kinematical 
variables for CP-odd observables: 

( be - Pe) • - ^) ) , 

(bt-A)-(^-^*)), (4) 

( [{Pe - Pe) X (Pt - Pt)] ■ - St) ) , 

where the spins and momenta are defined in the cm. frame. (The initial state is CP-even if we 
assume the SM interactions of with 7 and Z.) Generally one may think of other combinations 
involving Sg and as well. However, the spin directions of are not independent of their 
momentum directions for longitudinally polarized or unpolarized beams. Therefore, we would 
like to measure the difference of the spins (or the polarization vectors) of t and i. Practically 
we can measure the t and i polarization vectors efficiently using angular distributions. It is 
known that the angular distribution of the charged lepton from the decay of top quark is 
maximally sensitive to the top quark polarization vector. In the rest frame of top quark, the 
angular distribution is given by [^] 

1 dT{t M+u) _ 1 + Pcos^£+ . 
ft d cos 61+ ~ 2 ^ 

at tree level, where P is the top quark polarization and 6i+ is the angle of measured from the 
direction of the top quark polarization vector.[] Furthermore, we may think of CP-odd observ- 
ables bilinear in St and St, which require more complicated analyses for their reconstructions 
from decay products. 

We may anticipate following aspects of the CP-odd quantity 6P = (P — P)/2, (a half of) the 
difference of the t and t polarization vectors. It will be directly proportional to the coupling 
dt'y or dtz or dtg when only one of the couplings is turned on at a time. 6P will include a 
suppression factor (3 ~ \p^\/mt, the top quark velocity, since these couplings are accompanied 
by the top quark momentum; cf. eqs. @ and (^. Thus, this factor will be larger at higher 
cm. energy. Apart from this (3, energy dependence of SP originating from the anomalous tg 
coupling will be different from that originating from the anomalous ^7 and tZ couplings. The 
contribution of the coupling dtg to 6P will be suppressed if the energies of t and i are too 
large to allow for enhancement of QCD due to Coulomb binding effects. This happens if the 
cm. energy measured from the threshold £ = y/s — 2mt is much larger than the Coulomb 
binding energy, S ^ a'^rrit. On the other hand, the contribution of df^ or dtz coupling would 
not have such energy dependence since CP violation occurs at the first stage of the top pair 
production. Dependences of 6P on the e~ longitudinal polarization will be different between 
the photon-induced effect and the Z-induced effect, since ei and cr couple differently to 7 and 
Z. These differences in the energy and e~ polarization dependences can be used to disentangle 
the effects of the three anomalous CP violating interactions in the tt threshold region.^ 

II Indeed the distribution is ideal for extracting CP-violation in the ti production process; the above 
angular distribution is unchanged even if anomalous interactions are included in the tbW decay vertex, up to 
the terms linear in the decay anomalous couplings and within the approximation mt, — |Q . 

** Use of the longitudinal polarization for decomposing the photon-induced effect and the Z-induced effect 
was advocated in [po| . 
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3 The Top Production Vertices 



In this section we present the ^7 and tiZ vertices when the QCD binding effects and the CP- 
violating anomalous couplings are included. At tree level (and without anomalous interactions), 
the electroweak ti vertices are given by 



(r-)^ 



V 



A 1 



7 



Y75 (X = 7,^), 



(6) 



times —igx- Since the vertex (F"^)^ is contracted with the wave functions of 7 and Z produced 
by e'^e~ annihilation, only the space components of (T^)'^ are relevant. Here and hereafter, the 
Latin indices refer to the space components. See Appendix for the definition of the electroweak 
couplings Qx, v^^ and a*^. These vertices are modified by the QCD binding effects and by the 
anomalous interactions as 



(r^) =t;*^rv-a*^r^ + rf,xr 



tXl- X-EDM 



(7) 



where 



TT 



J m 



X 



Pi 

nit 



~{E + iTt 



X-EDM 



TT 



-ij,^F{E,p) xlP--{E + trt 



(8) 



p = Pt = ~Pt denotes the top quark momentum in the cm. frame, and p = \p\. We work in the 
potential-subtracted- mass scheme [^] instead of the pole-mass scheme, and E = y^— 2mps(/i/) 
represents the cm. energy measured from the twice of the potential-subtracted mass of top 
quark; rrit denotes the pole mass of top quark and it is expressed in terms of the potential- 
subtracted mass by 



mpsifJ^f) H /U/ 



TT 



1 + 



A-n 



ai-/3o(log^-2 



with 



ai 



(9) 



(10) 



where fif and /i denote the renormalization scale of the potential-subtracted mass and the MS 
coupling, respectively; Cp = 4/3, Ca = 3, Tp = 1/2 are the color factors and n/ = 5 is the 
number of active flavors. 

G{E,p) and F{E,p) are the S'-wave and P-wave Green functions, respectively, defined by 



— + Vir■,^^f)-iE + ^Tt) 
rrit 

^ + V{r-iif)-{E + iVt) 
nit 



G{E,x) = 5^(0;) , 



(11) 
(12) 
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and 



G{E,p) = j d^xe~'P '=G{E,x) 



(13) 
(14) 



V{r]fif) is the Fourier transform of the two-loop renormahzation-group-improved QCD po- 
tential, where the infrared renormalon-pole is subtracted and absorbed into the definition of 
mps(/i/). See for details. One may also write conveniently as 



G{E,p) = {p\ 



\x = 0) 



pF{E,p) = {p\—. —r^——=r-—wTP\x = 0) , 

p^/rrit + V — [E + zTt) 

The Green function associated with the gluon anomalous coupling is given by 

D{E,p) = G{E,p)-F{E,p). 



(15) 
(16) 

(17) 



We can see from eqs. (§) that the effects of all the anomalous CP-violating interactions 
are suppressed by \p\/mt ~ /?. Thus, for consistency we have incorporated all 0{as) = 0{P) 
corrections in the SM vertices. 

A sketch of derivations of these vertices goes as follows. Using non-relativistic forms of the 
t and i propagators. 



1+7° 



2 S/2 + ko -\k\y{2mt) +iTt/2 



1-7° 



2 Sl2-k^ -\k\'^/{2mt)+iTt/2 



(18) 



where = {2mt + S, 0) and k'^ = {k'^, k) in the cm. frame, a self-consistent equation for the 
vector vertex, similar to that given in p3, reads 



ryi£,pn = Y + CFi-tgsf 



d^A; 



(27r)4f/2 + P - \k\y{2mt)+iTt/2 
i 



X 



X 



8/2- k^- \k\^/{2mt) +iTt/2 



7 



-ry{8,k^)- 



7 



d-. 



2mt 



7°^^rV(^,A;)- 



X 



\p — k\ 



(19) 
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Since there is no p°-dependence on the right-hand-side, consistency requires ry{S,p'^) 
ry{S,p). Thus we can trivially integrate over and obtain 



7 



(27r)3 \p - fc|2 \k\ymt -{S + iVt) 



X 



X 



1 + 7° 



d 



1-7^ 



(20) 



We decompose the vertex function ry{S,p) into different spinor structures as 



1 + 0-° • 1 - 0-° 

^^ry{s,p)^-L 



7 



7TG(^,p)+7^(f^-^5^^)rB(^,p) 



IP 



1 



(21) 



By pluging this expression into the integral equation above, one obtains integral equations for 
scalar functions Tc{S,p), etc. One can see that Td{S,p) = 0{dtg), Tb{S,p) = 0(d^^^. Thus 
we neglect Vb{£,p) hereafter. Let us write 



Va{£,p)=[^-{£ + iT,)]G{£,p) 
Tp{£,p) = - {£ + iV,)] F{£,p) 

.2 



(22) 



Td{£,p) 



PL 



-i£ + iTt)]dtgD{£,p) 



Then G, F and D satisfy 



{£ + zT,)]G{£,p) 



El 

mt 

^_(£: + ,r,) )p^F(£:,p) + 

mt 

?^-{£ + iTt)]p^D{£,p) + 
nit 



d?k 
(2^ 
d^k 



{2nf 



Vc{\p-k\)G{£,k) 



Vc{\p-k\)k^F{£,k)=p^ 
Vci\p-k\)k'D{£,k) 



(23) 
(24) 



(27r)3 
d^k ~ 

^^^Vc{\p-k\)ik-pyG{£,k) 



(25) 



with Vc(g) = — Ci747ras/g2. Comparing the third equation with the first two equations, we find 
that D = G — F. The first two equations are equivalent to eqs. (|Tl|)- ([l6|) apart from the fact 
that the Coulomb potential —Gpag/r is replaced byQ V{r;fj,f) and that the renormalization- 
group-improved potential-subtracted-mass scheme |^ is used instead of the pole-mass scheme. 



* The replacement is justified: In Coulomb gauge the 0{as) corrections to the potential come solely from 
the vacuum polarization of Coulomb-gluon |40|. Hence, the net effect is to replace the fixed-coupling constant 
as(/i) in the leading-order by the V-schcme running coupling constant — k\). 
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The axial- vector vertex was derived in [p^j: 



mt 



(26) 



The hard- vertex factor for the vector vertex was derived in 411 and that for the axial- vector 



vertex in |^2|. We may also derive T^x-edm ^ similar manner. 

Two comments would be in order here. One might think that including the non-renormalizable 
interactions eq. (|l]) into loop integrals [e.g. eq. ([19|)] causes ultra-violet divergences and leads to 
unpredictability. We note that only non-relativistic domains of the loop integrals are relevant 
in resummations of the Coulomb singularities. In fact high momentum regions are effectively 
cut off in the self-consistent equations due to our non-relativistic approximation. Thus, we 
can calculate unambiguously the leading contributions of these effective interactions. In this 
regard, in eqs. (P) the hard-vertex correction factors are associated only with the SM contri- 
butions. We cannot determine hard-vertex corrections to the vertices including the anomalous 
interactions since non-renormalizability of these interactions matters at this order. 

The simple form of the Green function including the tg anomalous interaction eq. (p!7[) is a 
consequence of the following fact. The effect of Vcp-odd integrated over the time interval from 
t = to t = T can be written as {dtg/nit) {st — st) ■ [p{0) — p(T)] using the equation of motion 
p = — Wc. Namely the difference of the top quark momenta at t = and at t = T carries 
the net effect of the chromo-electric field which aligns the EDMs during < t < T. Concisely, 
for H = p"^ /rrit + Vq and 5H = Vcp.odd, the variation of the time evolution of the tt system is 
expressed as 

^t) ■ P 



mt 



-iHT 



Thus, the propagation at a fixed energy is given by {idtg/rrit) {st 
and ([1|). 



p{G 



(27) 

F); cf. eqs. (|TBp 



4 The Polarization Vectors and the Spin Correlation 
Tensor of t and t 



Using the vertices derived in the previous section we may write down the production cross 
section of a tt pair in the threshold region. The cross section, where {t,t) have momenta 
{Pt-i—Pt) and the spins +| along the quantization axes {st,St) in the cm. frame, is given by 

da{st,st) _ da 1 + P-st + P-St + {st)i{st)jClij 



d^Pt 



4 



(28) 



Here, \st\ = \st\ = 1. On the right-hand-side, da/d^p^ represents the production cross section 
when the spins of t and i are summed over: 

da / - ^ 



d^Pt 



1 



ACFas\ Nca'Tt 1-F,+ F,- 



TT 



X \G{E,pt)\\a, + xa2){l + 2Re 



''G{E,pt) 



(3 cos 6, 



te 



(29) 
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where /? = \p^\lmt and cos 6'te = PePtK\Pe\ ^ is the fine structure constant; Nc = ?> is the 
number of colors; x is a function of the initial longitudinal polarizations P(.±'- 



X 



1 - Pe+ Pe- 



rn 



If the positron beam is unpolarized (Pe+ = 0), x = ~Pe-- The coefficient Cfb and the constants 
oi, 02 are defined below. 

In eq. (|28|) P and P represent the polarization vectors of t and t, respectively. Both the 
SM and the anomalous interactions contribute to the polarizations: 



P = PSM + 5P 



PSM + 5P • 



(31) 



The SM contributions are CP-even (except for tiny CP-violating effects which we neglect) and 
are equal for t and i. On the other hand, the anomalous CP-odd contributions are opposite in 



sign: 



SM 



SM 



5P 



-5V . 



(32) 



Note that we are working up to linear terms in the anomalous couplings. Hereafter we express 
these vectors by components: 



P//n// + Plh^ + PnHn , St = s//n// + s±n± + snHn , 



(33) 



where the orthonormal basis is defined from the e beam direction and the top momentum 
direction: 



n// = 



Pe 

\Pe 



Pe^Pt 
\Pe^Pt 



n± = HNXn// 



(34) 



Then the polarization of t is given by 

(Psm)// = C° + Re (c}^^ PcosOte , 



(PsM)± = Re( Cx- )/3sin^te , 



SM JN 



Im ( Cn— ) P sin 9te 



for the contributions from the SM interaction^, and 

6P// = , 
5P, 



D 



Im ( BUtn^ 1 + Im ( BUt^^- ) + Im ( B^dtz 



5Pn 



D 



G 



F 
G 
F 



Re ( Bf,dtg- 1 + Re ( B^dt^- ) + Re ( B^dtz- 



G 



G 



f3 sin 6te 



(35) 
(36) 
(37) 

(38) 
(39) 

(40) 



^ These results were derived in 
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for the contributions from the anomalous interactions. The coefficients (7y etc. are defined 
below. 

There is also a term bilinear in St and St, which represents the correlation of t and t spins: 



(41) 



where 



and 



[st)iist)jQijgj^ = S//S// + {s//s_L + s±s//) Re ( Cn^ ) l3 sinOte 



+ (s//Sn + snS//) Im (^^^^^ 



(42) 



{st)i{st)j SQij = 

+ (s//Sn - SnS//) 



Im ( Bf^dtg^^ + Im (^B^dt^^^ + Im (^B^dtz^ 
Re ( BUtg^ 1 + Re ( B'^'dt^^ 1 + Re ( Bfdtz^- 



G 



G 



G 



(3 sin 6te 

PsinOte . (43) 



The coefficients Cj's and Bf^s included in Psm, ^P, QijsM' defined as follows. For 

the SM contributions,! 



c±{x) 



0-2 + XO-i 

ai + xa2 
04 + xas 



CUx) = 2C?,C. + C, = 2(1 - x Y;"\~"'':t , 



[ai + XO'2] 



<-'NlXJ — : — ltb 



The coefficients for the CP-odd contributions are given by 



(44) 



Bl{x)dt, + Bl{x)dtz 
B^ix) dt^ + Biix) dtz 



Q5 + XQ6 

0-1 + XO'2 

Q6 + XQ5 

Ol + X(^2 



(45) 



■I- Our notations are similar to those of |3l|,Q. There are two differences: (i) our 03 and 04 are factor two 
smaller than theirs; (ii) our Cat is defined in opposite sign to theirs. 
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and 

Blix) = -1 , 



ai + xa2 

Bfix) = { {[v'vr v^' + [aV]* a^^) + x {[v%r + K^T v^') } d{s) , 

ai + XO'2 

BUX) = CSix) , (46) 
ai + XO'2 



+ XO'2 

B^ix) = — ^ {X {[v%rv'' + [aV]*a^^) + ([t;V]*a^^ + [a%'Y v'')} d{s) . 

Symbols ai^g denote combinations of the electroweak parameters: 

ai = I [v%'] I ' + I [a'^t;*] \\ aa = 2 Re ( * [a^t;*] ) , (47) 

and 

= t;^^c/i^ + [v%'Y v'^d{s) dtz + [a^v'Y a^^d{s) dtz , 

ag = [v^vf [(^^d'] + [a%'Y l^^d'] (48) 

itz . 



fv'Y a^^diys) dtz + [a^v'Y v'^dt^ + [a'v'Y v'^d{s) dt 



Symbols [v'^a^] etc. are defined in the Appendix, eq. (^), and d{s) is the ratio of the Z 
propagator to the 7 propagator [eq. (|70D ]. 

We comment here why there is no ny component in 6P [eq. (^)] or why there are only 
a few components in 6Q,ij [eq. (^)]. SP and SQij originate from interferences of the leading 
SM amplitude Mo and the amplitude proportional to the anomalous couplings 6M. The SM 
amplitude Mq is in a linear combination of spin S// = ±1 states (c+ |TT)+C-|ii)), whereas 
the CP-reversed amplitude 6M is in spin S// = state (|Ti)~|ii))) see Sec. ^ In order to 
produce a non-zero interference between the two amplitudes, either one of the spins of t and t 
must be flipped. This is possible only by sandwitching the spin operator S^*^ or S^'' {i = t or 
t). 

We can understand from symmetry considerations the combinations of the electroweak 
couplings and the Green functions in each term of the production cross section da{st, St)/d^Pf.. 
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1 k 



-1 

-2 
-3 



B, 



perp Y 
perp-Z 



norm-g 
nonn-Y 



B, 



norm-Z 



-1 -0.5 0.5 

X % 

(a) (b) 

Figure 2: The electroweak coefficients Q's and i^f^'s (for P, P 
and Qij) vs. the initial polarization parameter x- In the figures, 



C, 



para 



C//, c, 



perp 



Cn, etc. (a) The coefficients for the 



SM contributions, (b) The coefficients for the contributions from the 
anomalous couplings. 



This provides a non-trivial cross check of the formulas presented in this section; the argument 
is presented in [H3 . 



5 Numerical Analyses of and SQij 

In this section we study numerically the polarization vectors and the spin correlation tensor 
derived in the previous section. We use the input parameters: mps(/i/) = 175 GeV, fif = 
3 GeV, /i = 20 GeV, mz = 91.19 GeV, a,(mz) = 0.118, and sin'^ 9w = 0.2312. 

First we examine the coefficients Cj's and Bf^s, which represent combinations of electroweak 
couplings. They are given as a function of x; cf. eq. (^0]). Fig. 0(a) shows the coefficients for 
the SM contributions Psm and Qijg^^- Except for C7y, typical sizes of the coefficients are order 
one. Fig. |^(b) shows the coefficients for the CP- violating contributions 6P and SQij. Typical 
sizes of all these coefficients are order one. We see that their dependences on x are different. 

Next we examine the S'-wave and P-wave Green functions. In Figs, ^(a) and (b) are 
shown the Green functions at -E = — 2 GeV and E = +2 GeV, respectively. They depend 
on both the energy E and the top momentum pt. The momentum distribution of top quark 
da/d^Pf. oc pj\G{E,pt)\'^ has a peak {pt = Ppcak) at a given cm. energy |2^, see Figs. ^ and|^. 



Then we may plot the ratios [3F/G and [3D/G, included in (5P and 5Qij, as a function of the 
energy E alone by choosing the top momentum to be the peak momentum Ppeak- These are 
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E = -2 GeV ; RG-improved NLO ; (i^ = 3GeV 



E = 2 GeV ; RG-improved NLO ; lif= 3GeV 




(a) 



1.5 
1 

0.5 


-0.5 
-1 
-1.5 





Re[G] 




Im[Gl 




Re[Fl - 




III. / , 






1 \J 



(b) 



Figure 3: The S'-wave and P-wave Green functions vs. the top quark 
momentum p at (a) E = -2 GeV, (b) E = +2 GeV. 
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Figure 4: p'^\G{E,p)\^ vs. the top quark momentum p at fixed cm. 
energies. These are proportional to the leading-order momentum dis- 
tributions of top quark. 




RG-improved NLO ; |J.y= 3GeV 




-4 -2 2 4 

E [GeV] 

Figure 5: The peak momentum Ppcak of the momentum distribution 
da/dp oc |pp|Gp vs. the cm. energy measured from the twice of the 
potential-subtracted mass, E = ^/s — 2mps{l^f)- It represents the 
typical momentum of top quark as a function of E. 
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At /) = /jpjjj; ; RG-improved NLO ; \lj = 3GeV 
0.2 I . • . . 




Figure 6: The ratios of the Green functions times the "velocity" of 
top quark evaluated at the peak momentum Ppeak of the momentum 
distribution da /dp; see Fig. ^. (a) These are given as a function of E. 
(b) These are plotted on a complex plane as E is varied. 



shown in Figs. |. We see that the size of \(3F/G\ is 5-20% while the size of \f3D/G\ is 5-10%. 
Clearly their energy dependences are different. Also it can be seen that the strong phases are 
quite sizable. 

One may understand these behaviors of the Green functions semi-quantitatively using an- 
alytic formulas. The relation 



Ppeak 



'mt {E + 1 GeV + iVt 



(49) 



agrees qualitatively with Fig. ^ in particular aX E > 0.| Here, 1 GeV ~ 2mt — Mis = 
"binding energy" . For a stable quark pair with the Coulomb potential, G and F can be 



obtained analytically for on-shell kinematics ||44||: 



lim 


(- 


^ 2 
Pt_ 






mt 








lim 


(- 


Pi 






mt 









+ iTt)G{E,pt: 



exp 



V=Vc 



V=Vc 



npB 
2pt 



T l + i 



■Pb 



1 — t — exp 

Pt 



Pt 

npB 
2pt 



r(i + .^ 

Pt 



(50) 
(51) 



where ps = CFOismt/2 ^ 20 GeV. Thus, we may find a sensible approximation formula 

Pb 



F 
G 



1-i 



P=PpGak 



^mt {E + 1 GeV + iVt) 



(52) 



This agrees qualitatively well with Figs. It follows that D/G ^ and F/G —>■ 1 when 
lE + iVtl > a^mt. 

Combining the analyses of the electroweak coefficients and the Green functions, we find 
that the typical sizes of the CP-odd quantities SP and SQij are 5-20% times the couplings 



For F ^ 0, Ppeak = ^/mt (E + iTt) 



holds exactly. 
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{dt^, dtz, dtg) in the threshold region. Using the different dependences on the polarizations 
and on the cm. energy, we will be able to disentangle the effects of the three anomalous 
couplings in the ti threshold region. A more comprehensive numerical study of the coefficients 
and the Green functions is presented in PSI . 



6 Observables and Sensitivity Estimates 

So far we have considered the CP-odd quantities related to the spins oft and t. These quantities 
are, however, not directly measurable observables in experiments. In this section we focus on 
5P and consider how to extract it. Then we give rough estimates of sensitivities to the CP- 
violating couphngs dt^, dtz, dtg expected at future e+e~ colliders. 

The top quark polarization vector can be extracted most efficiently using the angular distri- 
bution of charged leptons from the decay of top quarks. The charged lepton angular distribution 
in the ti threshold region, in the leading order, is given by 

da{e+e-^tt^U+iA)W-) ^ da 1 dTt^M+v{^) 

d^PtdVLi ~d3^^r\ dfi^ ■ ^ ' 

The left-hand-side shows that this is the differential cross section where the three-momentum 
of parent top quark and the direction of charged lepton are fixed, while all other variables are 
integrated over. The right-hand-side shows that it is given as a product of the ti production 
cross section eq. (|29|) and the decay angular distribution from free polarized top quarks. The 



top polarization vector is given by eq. (0). The lepton angular distribution in the top rest 
frame is given by eq. @. It coincides with the angular distribution in the laboratory frame 
in the leading order, since top quarks are almost at rest in the threshold region. Hence, the 
expectation value of the lepton three-momentum projected onto an arbitrary chosen direction 



n is proportional to the top quark polarization vector in the same approximation |31 



1 _l_ 2r + 3r^ 

((n-p^)) ^ 12(1 + 2r) ^ " " ^ ' (5^) 

where r = m^/m^, and ((■■■)) denotes an average taken for a fixed top three- momentum p^. 

Taking a CP-odd combination, the contributions of the anomalous interactions can be 
extracted as 

((n ■ (p^ + p^))) - ^ Q^^j^^^ mt X n ■ . (55) 



By choosing n = and Hn, we can extract the components eqs. (|39|) and (|40|) of 5P. The above 
formula remains valid even if we include the full 0{as) corrections (in particular the final-state 
interactions) in the SM parts of eqs. (^) and (|5^ ) |]3T|| , since the pure SM contributions 



drop in the CP-odd combination. Alternatively, we may consider a slightly different observable 

((n- (n, + n,))) ^ ^n-5P , (56) 

where n^/n^ denote the directions of This observable would be useful if the £^ directions 
can be measured more accurately than their three- momenta, such as in the case of r^. 
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In experiments the lepton-plus-4-jet mode can be used to reconstruct the lepton momen- 



tum and the top quark three-momentum simultaneously |2^. In order to detect a signal of 
CP-violation, it is not necessary to reconstruct the top quark three-momentum with a high 
accuracy. One should define the top quark momentum merely as the sum of all the visible 
momenta in a top-jet cluster, so no stringent cuts are required to reduce missing momentum. 
The only important point is that any experimental cut should be imposed in a CP symmetric 
way. Later when we measure accurately the values of the couplings dt^, dtz, dtg, we would need 
to reconstruct the top quark three-momentum to a reasonable accuracy. 

In order to extract 6Q,ij we need to measure spin correlations of t and i. Instead of eq. (|53[) 
we should consider a double differential decay distribution of t and i, which can be obtained 
using the formula of We may think of observables such as {{Pt — pt) ' (Pi ^ Pe)) ^^t CP-odd 
observables sensitive to 5Qjj. Here, we do not discuss extraction of SQij any further and leave 
the subject to our future work. 

Let us make rough estimates of sensitivities to the CP- violating couplings dt^, dtz, dtg 
expected in future experiments. Eq. ( ^Bf ) shows that a statistical reconstruction of the top 
quark polarization using lepton directions is quite efficient. The top quark polarization vector 
projected to a certain direction P = (s^ ■ n) is given by 

where {^i) denotes the number of top quarks with spin in the direction n (— n). Hence, 
the statistical error of P may be estimated as (5('^*'^*)p ~ l/^/Nes, where iVefj stands for the 
number of events used for the analysis. Assuming an integrated luminosity JC = 50 fb^^ and 
a detection efficiency e = 0.6, 

Nes = cTtt X JCx {2B,Bh) X e 

~ 0.5 pb X 50fb"^ X ■ ^ • x 0.6 = 4 x 10^ events, (58) 

which means 1/ ^ L5 x 10^^. The leptonic (hadronic) branching fraction i?^ {Bh) of 
into and //^ (hadrons) is given by 2/9 (2/3). 
Using the relations 



5P^ = Im ( Bldtg^ /3sin^te , ^Pn = Re (^^c?*,^ ) /^sin^^, , (59) 



the statistical error of dtg is estimated to be 



1 

d(cos6'te) 1 



1 



d(cos6'te) sm6'te 



1 4 
~ — X — X 1.5 X 10 

0.1 TT 



-1 

2 



0.2- (9(10%) . (60) 
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Note that we have similar sensitivities to both the real part and imaginary part of dtg inde- 
pendently using the two components of the top quark polarization vector. The above value 
translates to a sensitivity to the chromo-EDM of top quark at 

5^^*^*) (^rft,) r-W'gsCm. (61) 

Since all the electroweak coefficients {C±, i?N, etc.) are of similar sizes, and so are the ratios 
of the Green functions (Im(F/G), etc.), sensitivities to the EDM and Z-EDM are estimated 
to be at the same order: 



^(.tat) ( ±^ \ _ lO-i^ecm , ( ^d^) ~ IQ-^'g, 

\mt J \mt J 



cm . (62) 



A Monte Carlo simulation study is also in progress incorporating realistic experimental 
conditions expected at a future e~^e~ collider ||4^. They show that high detection efficiency is 



possible with simple event selection criteria and 6-tagging. Up to now, only the lepton energy 
asymmetry {Ei+ — E^-) for the dilepton-plus-2-jet events was studied The la statistical 
error corresponding to 100 fb~^ was obtained as 

They studied the bounds on the anomalous couplings setting the input values to be dt^ = 
dtz = dtg = 0. Based on our calculations, they obtained |Re[e*''^i'(ij^]| < 1.5, |Re[e*'^^(itz] | < 1.0, 
I Re [e*'^'' (if I < 3.9 at 95% confidence-level (statistical errors only), where e*'^-^'s are the relevant 
strong phases. In fact, the lepton energy asymmetry is not a good observable for extracting 
SP in the threshold region. It is suppressed by /5 ~ 10% compared to the CP-odd combination 
of the lepton three-memontum, eq. (^5]). Moreover the branching fraction for the lepton-plus- 
4-jet mode is larger than that for the dilepton-plus-2-jet mode. Thus, we expect that the 
sensitivities to the anomalous EDMs will be better by a factor 10 or more if we use the lepton 
three-momentum or the lepton direction. This is consistent with the naive estimates we made 
above. 



7 Summary and Conclusions 

In this paper we studied how to probe the anomalous CP-violating couplings of top quark 
with 7, Z and g in the tt threshold region at future e~^e~ colliders. The anomalous couplings 
contribute to the difference of the t and i polarization vectors, 6P = (P — P)/2, as well as to 
the spin correlation tensor 6Qij. We studied dependences of these CP-odd quantities on the 
beam polarizations, cm. energy, and top quark momentum. We find that the typical sizes of 6P 
and 6Qij are 5-20% times the couplings {dt^,dfz,dtg) in the threshold region. Experimentally 
we can measure SP efficiently using the expectation value of the CP-odd combination of the 
£^ momenta, Pi + P£, or of the £^ directions, -|- n^. We have similar sensitivities to both the 
real part and imaginary part of dt-y, diz, dtg independently using the two components of the 
top quark polarization vector 6P± and (5Pn. Taking advantage of different dependences of SP 
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on the polarizations and on the cm. energy, we will be able to disentangle the effects of the 
three anomalous couplings (ij^, d^z, dtg in the ti threshold region. We made rough estimates 
of sensitivities to the anomalous couplings expected at future e"'"e~ colliders, considering as a 
simplest example extraction of 5P from the £^ distributions. For an integrated luminosity of 
50 fb~^, we estimated 

5^''^'^dtz, 5^''^'^dtg ~ 0(10%) (64) 

when only one of the couplings is turned on at a time.|^ The above values translate to sensi- 
tivities to the top quark EDM, Z-EDM and chromo-EDM: 



^(stat) I ^dtA ~ 10-i^ecm , (J^*^*^*) \ —dtz ) ~ 10 ^'fecm 
mt J \mt 



^(stat) / i^dtA ~ 10^1^^, cm . 
\mt J 



(65) 



The sensitivities to the top quark EDM and Z-EDM are comparable to those attainable in 
the open-top region at e^e~ colliders |T^. The sensitivity to dtg is worse than that expected 



at a hadron collider P,^p!0|,p!3|-[16|| but exceeds the sensitivity in the open-top region at e^e~ 
colliders We note that there is an advantage of the ti threshold region. The clean 

environment of an e^e~ collider enables accurate determination of the value of the top-gluon 
anomalous coupling dtg if its value happens to be larger than 0(10%). On the other hand, 
at hadron colliders it would be difficult to measure the value of the coupling with a similar 
accuracy even if a CP-violating effect is detected. Regarding energy upgrading scenario of a 
future linear e^e~ collider, it is possible that the machine operates first in the ti threshold 
region for a significant amount of time, while measuring the top quark mass precisely, etc., 
before the beam energies will be increased to the open-top region. Therefore it would be 
desirable that measurements of the anomalous couplings can be performed concurrently with 
other unique measurements near threshold, with sensitivities comparable to those in the open- 
top region. Unfortunately the sensitivities to the CP-violating couplings achievable in the ti 
threshold region are one to three orders of magnitude larger than the predicted sizes of top 
quark EDMs in the models reviewed in Sec. |l|. Using the results of this work, a Monte Carlo 
study incorporating realistic experimental conditions expected at a future e^e~ linear collider 



is underway 46 1. 
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Appendix: Conventions and Notations 



In e^e tt, both 7 and Z are exchanged in s-channel. Their effects can be combined in 
terms of effective couplings. We denote the SM vertices for electron and top quark by 

Ax^ = gx {v'^'l, - a^^T^Ts) , = 5,, {v'^'r - a'^'n^) {X = 7, ^) , (66) 

where 

— e — gsinOw , v^'^ — Qf , a^'^ — , 

9z = ^^, v^^ = Itsl- Qf sin' 9w , a^^ = ]-T,l ■ (67) 
cos6'h^ 2 2 

The amphtude for e^e" — > tt at tree level can be written as 

o Tit V" 

X=y,Z ^ 

e2 r 

= — b'f^*] {'"iPehnUipe)) (uipth^vipt)) 
S L 

- [V^a^] {v(pe)-fij,u(pe)) {u{pt)l''lbV{pt)) 

- [(^"A {v{Pe)lnlbU{Pe)){u{pt)l''v{pt)) 

+ [a^a*] {v{pe)-f;,j5u{pe)) {u{pth^j5v{pt)) , (68) 

where 



[v%'] 


= -y'^f;*^ 


+ d{s) 




1 






= v^^a*^ 


+ d{s) 




^d{s) 






= a'^v''' 


+ d{s) 




= d{s) 






= a^^a*^' 


+ d{s) 




= d{s) 





(69) 



represent energy- dependent "couplings" . Extensions to the anomalous vertices should be ob- 
vious. d{s) is the ratio of the Z propagator to the 7 propagator 

with 

9 



{g'm\Q-wY cos^^^ysin^^ 



= 5.625 . (71) 



w 

The width Y z of Z introduces an absorptive part. At a/s = 2 x 175 GeV, its relative magnitude 
is s/(s — vp?z + ^'^z^ z) = 1.073 — 0.002 i. Thus, in the threshold region, the Coulomb binding 
effects overwhelm the effect of Tz- Also in the open-top region, it is known that the QCD 
correction is larger than the effect of F^, as far as the normal component of the polarization 
of top quark Pn is concerned. 
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